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dj ! Abstract 
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^ ■ We apply a simple mean field like variational calculation to compact QED in 2+1 di- 

mensions. Our variational ansatz explicitly preserves compact gauge invariance of the 
theory We reproduce in this framework all the known results, including dynamical mass 
generation, Polyakov scaling and the nonzero string tension. It is hoped that this simple 
example can be a useful reference point for applying similar approximation techniques to 
nonabelian gauge theories. 
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1. Introduction. 



Study of the confining regime in QCD, as well as of many other strong interaction 
phenomena requires going beyond the perturbation theory. However, application of ana- 
lytic nonperturbative methods in quantum field theory, is a very complicated and not too 
well developed area. This is especially true for nonabelian gauge theories. Recently we 
have formulated a gauge invariant variational approach and used a restricted variational 
ansatz to study the ground state of a pure Yang Mills theory in 3+1 dimensions [|TJ. As 
with any new technique, it is desirable to develope an intuition for it, by first considering 
simpler examples. 

A theory which posesses many common qualitative features with QCD (such as confine- 
ment and dynamical mass generation), and yet is much simpler and much more tractable 
is compact electrodynamics in 2+1 dimensions. Moreover, this theory has been previous 
extensively studied by both, analytical || and numerical || methods. It seems therefore 
to be a perfect testground for application of our variational method. This is precisely the 
aim of this note. We will apply a gauge invariant variational approximation of ref.JT] to 
this theory. It is hoped that this toy calculation can teach us something about improving 
variational ansatz for realistic 3+1 dimensional nonabelian theories. It is also a nice excer- 
cise in itself, since it gives a vivid hamiltonian picture of Polyakov's monopole-instanton 
condensation phenomenon, which to our knowledge does not exist in the literature. 

The paper is organized as follows. In the rest of this section we discuss the Hamiltonian 
formalism for the compact QED3. In section 2 we set up our variational ansatz and 
discuss some important properties of the variational wave functionals. Section 3 contains 
calculation of the expectation value of the energy and solution of minimization equations. 
We also show that the Wilson loop in the best variational state has the area law and 

2 



calculate the string tension. In section 4 we discuss our results and the interpretation of 
our calculation from the point of view of Polyakov's dilute monopole gas approximation 
to Euclidean partition function. 

The theory is defined by the following Hamiltonian: 

H= 1 -[E^ + b 2 ] (1.1) 

The field b is somewhat different from the usual magnetic field B = eijdiAj. We will 
explain its definition in a short while. All the physical states should satisfy the Gauss' 
law constraint 

exp{i J d 2 xd i X(x)E i (x)}\m >= |tf > (1.2) 

One should note that there is a crucial difference between the Gauss' law in the compact 
theory and in the noncompact one. In the noncompact theory equation (|1.2|) should be 
satisfied only for regular functions A. For example, the operator 

which has the form of ( |1.2| ) with the function A proportional to the planar angle 9, 
A = ^9(x), does not act trivially on physical states. In fact, this operator creates pointlike 
magnetic vortices with magnetic flux 27r/g Q and therefore changes the physical state on 
which it is acting. 

In the compact theory the situation in this respect is quite different. Pointlike vortices 
with quntized magnetic flux 2nn/g can not be detected by any measurement. In Euclidean 
path integral formalism of ref. this is the statement that the Dirac string of the monopole 
is unobservable and does not cost any (Euclidean) energy. In the Hamiltonian formalism 
this translates into the requirement that the creation operator of a pointlike vortex must 
be indistinguishible from the unit operator. In other words, the operator ( |1.3| ) generates 
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a transformation which belongs to the compact gauge group, and should therefore act 
trivially on all physical states. Equation (|1.2|) should therefore be satisfied also for these 
operators. 

Accordingly, the Hamiltonian of the compact theory also must be invariant under 
these transformation. The magnetic field defined as B = eijdiAj, on the other hand does 
not commute with V(x), 

2-7T 

V\x)B(y)V(x) = B{y) + — 5 2 {x - y) (1.4) 

9 

The Hamiltonian should therefore contain not B 2 but rather only its singlet part. This 
is the meaning of the field b in equation (|1.1| ). Formally 

b 2 = PB 2 P (1.5) 

where P is the projection operator on the whole compact gauge group, which includes 
V(x). This form is convenient for the purposes of the present calculation, and we therefore 
do not write down a more explicit expression for 6 2 Q 



2. The variational ansatz. 



Our aim in this paper is to find a vacuum wave functional of this theory. Following 
Polyakov, we will be working in the weakly coupled regime. Since the coupling constant 
g 2 in 2+1 dimensions has dimension of mass, weak coupling means that the following 



*In the standard Euclidean lattice formulation of the theory the potential energy has the form 
cos(ga 2 _B), where a is the lattice spacing. This obviously has the same property as 6 2 , that is invariant 
under the transformation eq. 1.4 and reduces to it in the weak coupling limit. Our definition eq. 1.5 is 
equivalent to the Villain form of the action. 
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dimensionless ratio is small 

C « 1 (2- 1 ) 
A 

Here A is the ultraviolet cutoff in the momentum space, which as always has to be intro- 
duced to regularize a quantum field theoryfj. For a weakly coupled theory one expects the 
vacuum wave functional (VWF) to be not too different from the vacuum of a free theory. 
Since the VWF of free (noncompact) electrodynamics is gaussian in the field basis, 

*l>[Ai] = exp {-1/2 J #x#yAi{x)G£{x - y)A J {y)} (2.2) 

Gaussian variational approach in this case should give a good approximation. An impor- 
tant caveat, however is that the ground state WF should be gauge invariant under the 
full compact gauge group. As a result it turns out that one can not take just a Gaussian 
in A4, since this will not preserve gauge invariance. The simplest generalization of the 
Gaussian ansatz which we use along the lines of ffl], is to project a gaussian WF into the 
gauge invariant subspace of the Hilbert space. We therefore take as our variational anzatz 
the following set of states 

A i {y)--d i (j> i {y)) 
9 JJ 

(2.3) 

The functional integral is over the phase function 4>(x), and correspondingly the derivatives 
of (f) in the exponential are understood modulo 2n. That is, these derivatives do not feel 
quantized discontinuities in <p(x). The mathematically more precise way to write this 
is to substitute <9j0(x) by exp{— i(j)(xj}di exp{z</>(x)}. We will use however, the above 

shorthand notation for convenience. The simple rotational structure of Gij = 5ijG that 

*If one considers compact U(l) as an unbroken sector of the spontaneously broken SU(2) (Georgi- 
Glashow model), the ultraviolett cutoff A will be proportional to the scale of the symmetry breaking. To 
be more precise it is by order of magnitude of the charged vector boson mass Mw- 



^[Ai] = J L>0exp j-iy d 2 xd 2 y . 



Ai(x) - ~di<j)(x) 
9 



G-\x-y) 
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appears in the variational wave functional (|2.3|) is consistent with perturbation theory, as 
discussed in \M. 



The ansatz (|2.3|) depends on one function G(x). We now have to calculate the expec- 
tation value of the energy in this state, and then minimize it with respect to G(x). 

Before proceeding with the calculation we make the following (obvious) comment. The 
trial wave functional (|2.3| ) has a simple interpretation from the point of view of states of 
a noncompact theory. To see this let us rewrite the functional measure Dip in a slightly 
different way. Any angular function 4>(x) can be parametrized as 

<j>{x) = (j>(x) + 4>v{x) (2.4) 

where 4>(x) is a smooth function and <p v (x) contains all the discontinuities and can be 
written as 

Mx) = T,0(x-Xa)- y E i e ( X - X P) ( 2 ' 5 ) 

o=l (3=1 

and 9(x — x a ) is a polar angle on plane with a center at x a . The functional measure can 
be written as 

00 00 

D$ E E E E (2-6) 

n+=l n_=l {x a } {xp} 



00 00 1 « ■ iii— , . 

^EE -rn nn / fcafxpA' 

n+=0n_=0 "'+•"'-• a= ip = i V- 7 



where in the last equality we have just substituted integration over the coordinates of 
vortices and antivortices for summation, and for this reason introduced explicit UV cutoff 
A. We also introduced term with n + = n_ = corresponding to the absence of vortices. 

Let us define the function 

[ ' G-\x-y) 



Ai(x) - ^di<p(x) 



Mv) - -diMv) 
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(2.7) 



which differs from \1/ in that the integration is performed only over continuous gauge 
functions <fi. Obviously, ^ is invariant under noncompact gauge group, and therefore 
belongs to the Hilbert space of the noncompact theory. When acting on it, the vortex 
operator V(x) defined in (|1.3|) just shifts Ai(y) by ^d{6(x — y). We have therefore the 
following representation for *&[A] 

oo n + n_ 

*M = E UUVMV*(xp)9[A] (2.8) 

n+,n_=0 a=l /3=1 

This representation makes explicit the fact that a WF of the compact theory is constructed 
from a WF of the noncompact theory by taking a superposition of arbitrary number of 
vortices and anti vortices at every point. This superposition is obviously invariant under 
multiplication by a vortex operator and therefore is its eigenfunction with eigenvalue 1. 

Having noted this, we now proceed to calculation of expectation values in the trial 
state (|2.3|) and the minimization of the vacuum expectation value of energy. 



3. The energy minimization. 



Let us for convenience introduce the following notation 

At{x) = Mx) -~d4{x) 
9 

We will also switch to the matrix notations in the following, so that 



AiMAt = / d 2 xd 2 yA i (x)M(x-y)A i (y) 



(3-1) 



(3.2) 



The expectation value of any operator in the WF ( |2.3[ ) is calculated as 

< O >= Z- 1 J D(j)'D(j)"DAi exp {-^f G^A? } 0(A) exp {-^f G~ x Af j (3.3) 



Here Z is the normalization factor, which is just the norm of the trial wave functional 
(|2.3| ). Further, if the operator O is explicitly gauge invariant, we may shift the integration 

Aft 

variable A t — > A? , and reduce this expression to 

< O >= Z- 1 J D(f)Dr]DAi exp i~M G~ l Af X 0(A) exp {-^G^a} (3.4) 

where we have defined = 0' — <fi" and 77 = + <fi" . The integral over 77 is trivial and just 
gives the volume of the gauge group. Since the same integral exactly enters Z, it always 
cancells between the numerator and denominator, and we shall omit it in the following. 

As a first step let us calculate the normalization factor Z. 

Z = J D(j)DAi exp {-- [AfG^Af + A i G~ 1 A l ] j (3.5) 

The integral over A4 is gaussian and can be trivially performed. The integral over the 
noncompact part of the gauge group <fr is also gaussian. Then 



Z = Z a Z^Z v (3.6) 



with 



Z a = det(TrG) (3.7) 
= jDj) exp |-^2 / d^G-'d^ = det ' " 1 

Z v = J D(p v exp J di^G^diCpA 



Using equations fl2.5| ) and fl2.7|), one can represent Z v as a partition function of the 
$as of vortices 

00 n + n_ 

Z v = E 1111 dx a d Xf3 z n + +n - (3.8) 

n+,n_=0a=l (3=1 

exp < - j^E D ( x <* ~ Xa ') + H D ( X P ~ X P') -^ D (x a - xp)] 

[ 9 a,a' I3,j3' a,f3 



The vortex - vortex interaction potential D(x) and the vortex fugacity z are given by 

D(x) = 8tt 2 J 'J^LG-\k)cos(kx) (3.9) 

z = A 2 exp{-^D(0)} 

Here G(k) is the Fourier transform of the variational "propagator" G(x). 

The formula (|3.9|) reminds one of Polyakov's partition function of the monopole gas 
||. One should keep in mind, however that the physical meaning of it is quite different. 
The gas described by equation ( |3.9|) is two dimensional and not three dimensional, and the 
interaction between the particles is not Coulomb, but rather depends on the variational 
function G. Nevertheless, it has a definite relation with Polyakov's gas of monopoles and 
we will discuss this point in the last section of the paper. 

Since D(0) is singular, the last equation should be understood, as usual in the regular- 
ized sence, that is at finite UV cutoff, D(Q) should be substituted by D{x = 1/A). Note 
that the variational function G explicitly appears in the vortex - vortex potential. Since 
we expect the UV behavior of G(k) to be the same as in the free theory (G^ 1 (k) — > k), 
we have 

z = A 2 exp{---} 3.10 
2 g l 

In the following we will need to calculate correlation functions of the vortex density. 
To facilitate this we use the standard trick ||, || to rewrite the partition function Z v in 
terms of a path integral over a scalar field \- Let us introduce the vortex density p(x) 

p(x) = 8(x — x a ) — 8(x — xp) (3-11) 

The exponential factor in eq. (including the factors of fugacity) can then be rewritten 

as 

^2(71+ +n_) J Dxexp{ _2g 2 X D- 1 x + i px} (3.12) 
9 



The summation over the number of vortices is trivial and gives 

Z v = J D X exp{-2g 2 X D- 1 x + J^A 2 C os X (x)} (3.13) 

To calculate the correlator of p one can add ipJ to the vortex free energy, and calculate 
functional derivatives of the resulting partition function with respect to J at zero J. A 
simple derivation gives 

< p(x)p(y) >= Ag 2 D-\x - y) - lQg 4 < D- l x{x)D- l X {y) > (3.14) 

The propagator of x is easily calculated. At weak coupling z is very small, and all 
our calculations will be performed to first order in z. To this order the only contribution 
comes from the tadpole diagrams. This is easily seen by rewriting the cosine potential in 
equation ( |3.13| ) in the normal ordered form * 

2; 

cosx =< cosx > : cosx := ■ cosx ■ (3.15) 
Therefore to first order in z 

< *M*(0) 4g!D _! {k) + 2: = ^ _ z £^> + ^ ,3,6) 
The correlator of the vortex densities is then 

K(k) = f d 2 xe ikx < p{x)p{0) >= 2z + o{z 2 ) (3.17) 



and in this approximation does not depend on momentum, the fc-dependence will appear 
in z 2 and higher order terms. 

Now we are ready to calculate the expectation value of the Hamiltonian ( |1 . 1| ) . First, 
consider the electric part 

< / (l 2 .rE 2 > - < i d 2 x-^-2 > 



"The normal ordering is performed relative to the free theory defined by the quadratic part of the 



action in equation (3.13) 
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Z~ x I D^DAiexpl-^AfG^At 



2trG 



-i 



AiG~ 2 A 



exp 



1 



-AG -1 A 



(3.18) 



= trG- 1 - -^Z^Z- 1 J D0 (d^G- 2 ^) exp ^--Lo^G^d^ 
Performing the gaussian integral over 0, this reduces to 



(2 
d 2 k 



r 

.G-\k)-%z 
9 



(3.19) 



d 2 k 



k- 2 G- 2 (k) 



Let us note that it is only because of vortices one has negative contribution to the electric 
part of the energy. 

Now for the magnetic part. Since b 2 is the singlet part of B 2 , by definition in every 
gauge invariant state < b 2 >=< B 2 >. We will therefore calculate < B 2 >. Here one 
should be a little careful. Since B 2 itself is not gauge invariant, one can not use equation 
but rather explicitly keep both integrals, over and r\. 

1 



<b z > 



z- 1 
z- 1 



"DAB 2 exp 
D(j)Dr]DA 



AfG^Af + A\ G~ X A\ 

\ 1 2 



(3.20) 



AfG^Af + AG -1 A 



where, the factor Z contains an extra factor of the volume of the gauge group relative 
to equation (|3.5| ), and as previously = 0' — 0" and rj = 0' + 0" . The linear term in rj 
vanishes due to the symmetry of the measure under transformation rj — > —77. The term 
quadratic in rj is independent of the variational parameter G. It does not contribute to 
the minimization equations, and we omit it in the following. We therefore obtain 

/ \ -1 2 



< b 2 > = Z- 1 D(j)DrjDAi 



eijdi I Aj - — c 



exp 



1 



Z' 1 J DAiltijdiAj} 2 exp{-AG _1 A} = i J d 2 kk 2 G(k) 



(3.21) 



11 



This is the same result as in the noncompact theory. In fact, this is precisely what one 
expects, since a compact state \1/ differs from a noncompact one, ^ only by the presence 
of vortices, but b 2 by definition should not feel their presence. 

Summarizing, we have the following expression for the expectation value of the energy 
density 



1 1 r d 2 k 

v <h>= aJ (2^ 



4vr 2 



G-\k) + k 2 G(k) - —zk- 2 G- 2 {k) 



r 



The minimization equation is 



k 2 - G' 2 {k) 
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2zk- 2 G- 3 (k)-47i 2 -^- ( -fP- p - 2 G- 2 (p) 



5G(k) J (2tt 



(3.22) 



(3.23) 



From equation ( 3.10Q one finds 



5z 



5G(k) 



1 



-^k- 2 G- 2 (k)z 



Assuming perturbative behaviour of G at large momenta (G(k) 
the fourth term in equation ( 3.23|) to the third term is of order 



Sz Jd 2 pp- 2 G- 2 (p) A 2 
5G(k) 2zk- 2 G- 3 {k) X ~cfk 



(3.24) 



k 1 ), the ratio of 



(3.25) 



At weak coupling this is much greater than one for any value of momentum. We can 



therefore omit the third term from equation Q3.23|) and get a very simple equation on 

G- 2 (k) 



k 2 - G- 2 {k) = ^zk~ 2 G-\k) J ^T 2 P- 2 G~\P) 



with the solution 



where 



9' 



G- 2 (k) 



k 4 



k 2 + m 2 



4vr 4 r d 2 k 



m 



9 4 J (2vr) 2 
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k- 2 G- 2 (k) 



(3.26) 



(3.27) 



(3.28) 



Using equations ( |3.10|) and (|3.27| ) we get an equation on mass 



2 4tt 4 / vr 2 f d 2 V 1 \ f d 2 k k 2 o 

m = _a exp y , „, 2 i ) / ,„„, 2 ,„ 2 , ^ 2 (-^>) 



g 4 ^\ g 2 J {2tt) 2 y/p 2 + m 2 J J (2tt) 2 k 2 + m 2 
which in case of weak coupling A/g 2 >> 1 can be simplified as 



m 2 = vr 3 — exp{---} (3.30) 



,A 4 r ttA 

and because m « g 2 << A one indeed could neglect the m-dependence in the right 



hand side of equation ( |3.29|) . It is clear, that m is precisely the mass gap of the theory. 



Calculating, for example the propagator of magnetic field, we find 

J d 2 xe ikx < b(x)b(0) >= ^(k 2 + m 2 ) 1 / 2 (3.31) 

Note, that the dynamically generated mass we obtain in our approximation agrees 
with Polyakov's result @. 

Perhaps the most interesting question is whether our best variational VWF is confin- 
ing. To answer this question we calculate the expectation value of the Wilson loop. 

W c =< exp{ilg I A^} >=< exp{ilg f BdS} > (3.32) 

where I is an arbitrary integer and the integral is over the area S bounded by the loop 
C. We have written B rather than b, since this exponential operator is invariant under 
transformations B(x) — > B(x) + — , generated by the vortex operator. 



W, 



c 



Z' 1 J DAiexpl-AiG^Ai + ilg JetjdiAjcPxl (3.33) 
Z^Z- 1 J D0exp | -^di^G~ x d^ + z~ jf drfdxi 1 



The first factor in weak coupling is simply 



W = exp l-^9 2 J xy < B(x)B(y) > d 2 xd 2 y^ (3.34) 
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where the integral over both x and y is over the area S. In the limit of large S the leading 
piece in the exponential is 

^Jm^)^ 3 2 mS (3.35) 

This term therefore leads to the area law behaviour and gives the string tension 

I 2 

a = —g 2 m (3.36) 

and we see that area law Wo ~ exp(— aS) is a direct consequence of the nonzero mass 
gap m. 

The second factor in equation (|3.34j) is different from unity only for odd Z, since 
§ di<fidli = 27c(n + — n_), where n + (n_) is the number of vortices (antivortices) inside the 
loop. For odd I it can be easily calculated 

W v =< exp{iir J s p{x)d 2 x} >= J D X exp{-2g 2 xD' 1 X + J 2A 2 cos(x(x) -a(x))} (3.37) 

where a(x) is the function, which vanishes for x outside the loop, and is equal to tt 
for x inside the loop. At small coupling this can be calculated in the steepest descent 
approximation. The solution to the classical equations which contributes to the leading 
order result is x{ x ) = 0- For this solution 

W v = exp{-4^5} (3.38) 

where we again used the normal ordering prescription as in ( 3.15|) . Clearly this is a 



subleading correction to string tension ( [3.36D , since z oc m 2 and z/g 2 m ~ Jz~j 1 K 



exp(— irA/4g 2 ) << 1. With this exponential accuracy, the string tension is therefore 
given in our approximation by equation (|3.36j ). 
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4. Discussion. 



We have presented a simple variational calculation in the compact QED 3 . Our vari- 
ational ansatz was the direct adaptation of the ansatz of ref.fll] to this theory. The trial 
wave functionals are explicitly gauge invariant under the compact gauge group. The in- 
tegration over the gauge group is directly responsible for nontrivial dependence of the 
energy expectation value on the variational parameters, which leads to the generation of 
scale in the best variational state. The correlators and the Wilson loop calculated in the 
best variational state agree with known results. 

It is illuminating at this point to interpret our calculation from the point of view of 
the three dimensional Euclidean path integral. The vacuum wave functional of the theory 
can be represented in path integral formalism. To get the vacuum WF *&[A] one should 
calculate the path integral over the fields A(x,t), with t varying from — oo to 0, with 
the boundary condition A(x,t = 0) = A(x). To be more precise, in calculating VEV of 
some operator 0(t = 0), one should split the time coordinate of the plane with the time 
coordinate of the operator, so that one considers ^f[A(x,t = — e)] and Hf*[A(x,t = e)] in 
the limit e — > 0. 

The basic objects that appear in the Euclidean path integrals are monopoles, which in 
3D are not propagating particles, but rather instantons. When described in terms of the 
vector potential, or noncompact field strength, a monopole has a Dirac string attached 
to it. It is clear that the vortices (antivortices) of the gauge function 4>(x), in equation 
( |2.3| ) correspond precisely to the intersections of the Dirac strings of the 3D monopoles 
(antimonopoles) with the equal time plane at t — 0. The positions of the Dirac strings 
are not physical in the compact theory, and only the position of the monopole itself is 
gauge invariant. In fact, for all monopoles that do not sit in the infinitesimally thin time 
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slab between the planes t — — e and t = e, one can always choose the direction of the 
Dirac string such that it does not intersect the two planes. This precisely corresponds to 
expression equation (|3.4|) . The combination that enters this path integral nontrivially is 
= 0' — 0". At the points, where both functions 0'(x) (which corresponds to and 
cf)"(x) (which corresponds to have a vortex, 0(x) is regular. This is the situation, 
when a Dirac string intersects both planes t = ±e. When a 3D monopole sits in the slab, 
only one of the functions 0' or 0" has nonzero vorticity, and so does 0. The integration 
over 0(x) in equation ( |3.4| ) can be interpreted therefore as the direct contribution to the 
expectation value due to the monopoles at precisely the time t = 0. 

The fact that in this way one sees directly only the monopoles at t — 0, does not mean 
of course that other monopoles are not taken into account in this approximation. Indeed, 
the "bare" interaction potential between the t = monopoles is D(x) of equation ( |3.10| ). 
In the best variational state it is already short range, as follows from the solution for G(k) 
(see equation ( |3.27p ). This is in accordance with the 3D picture, that the 3D monopole 



gas produces screening. Obviously, if one only looks at the thin slab, every monopole there 
will have an antimonopole partner, which sits nearby (inside the screening length) in the 
third direction. The 2D monopole gas will therefore be screened by the 3D interaction, 
even before the interaction of the 2D monopoles between themselves is taken into account. 
This is perfectly consistent with our calculation. It is interesting to note, that even though 
this 2D interaction produces additional screening (the cosine term in the effective theory 
(|3.13|) ), it is the 3D screening that is responsible for the area law of the Wilson loop, as is 



clear from the calculation of the string tension in (|3.34j ). In fact, if one takes for G(k) the 



noncompact expression, both, the leading part Wo, equation (|3.34j) , and the subleading 



part W v , (|3.37|) have the perimeter law behaviour. The string tension in Wq vanishes, 



because in this case lim fc ^ k 2 G{k) = 0. In W v in this case D 1 {k 2 = 0) = 0, and for 

16 



large loops the classical equations of motion, which follow from equation (|3.37|) , apart 
from x — have another solution, which leads to the perimeter dependence of W v . The 
existence of the cosine term in the interaction (which appears due to the 2D screening) 
does not preclude the existence of this extra solution. 

Finally, we note that the calculation presented here can be extended to compact 
QED in 3+1 dimensions. In this case the vortex gas part of the partition function will 
be replaced by the gas of vortex loops with the interaction between the loop elements 
dependent on the variational function G. It is well known that compact QED^ has a 
phase transition at a finite coupling constant. There is a good chance that this phase 
transition will be seen in the present approximation. Consider variational functions G, 
for which d 2 G 2 is short range, like in equation (|3.27|) . The interaction between the strings 
is then short range. In this case the standard energy vs. entropy argument is telling us 
that at large g 2 the strings are condensed. The string gas contribution to the vacuum 
energy will then be sizable. Since this contribution has a negative sign, this situation will 
be energetically favored. Therefore at large g 2 one expects the best variational state to 
have a short ranged d 2 G 2 , and therefore the mass scale will be generated dynamically. At 
small g 2 the vortex rings do not condense even for short range interaction between them. 
Their contribution to the vacuum energy will be negligible, and one expects that the best 
variational vacuum will be determined by the contributions from the Ai and <fi integrals, 
which will lead to the same solution as in the noncompact theory. 
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